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ABSTRACT 

We present the first robust evidence of an anti-correlation between the X-ray photon 
index, T, and the X-ray luminosity in a single low luminosity active galactic nuclei 
(LLAGN), NGC 72f3. Up to now, such anti-correlation trends have been seen only in 
large samples of LLAGN that span a wide range of X-ray fluxes, although the opposite 
behaviour (i.e. a positive correlation between L and X-ray luminosity) has been ex- 
tensively studied for individual X-ray bright active galactic nuclei. For NGC 7213, we 
use the long-term X-ray monitoring data of Rossi X-ray Timing Explorer (RXTE), 
regularly obtained on average every two days from March 2006 to December 2009. 
Based on our X-ray data, we derive the T versus flux and the hardness ratio versus 
flux relations, indicating clearly that NGC 7213 follows a 'harder when brighter' spec- 
tral behaviour. Additionally, by analysing radio and optical data, and combining data 
from the literature, we form the most complete spectral energy distribution (SED) 
of the source across the electromagnetic spectrum yielding a bolometric luminosity of 
1.7 x 10 43 erg s _1 . Phenomenologically, the SED of NGC 7213 is similar to that of low- 
ionization nuclear emission-line region. The robust anti-correlation trend that we find 
between T and X-ray luminosity together with the low accretion rate of the source, 
0.14 per cent that of Eddington limit, make NGC 7213 the first LLAGN exhibiting a 
similar spectral behaviour with that of black hole X-ray binaries in 'hard state'. 

Key words: galaxies: individual: NGC 7213 - X-rays: galaxies - galaxies: nuclei - 
galaxies: Seyfert - X-rays: binaries - accretion, accretion discs 



1 INTRODUCTION 

NGC 7213 (z=0.005839) is a face on Sa galaxy hosting an 
active galactic n ucleus (AGN). It has been classified as a 
type I Seyfert bv lPhillipsl l| 19791 ). based on its Ha linewidth 
(full width at zero intensity of 13000 km s -1 ), but also as 
a low-ionization nuclear em ission-line region (LINER) by 
iFilippenko fc Halpernl (|l984l V based on the study of a vari- 
ety of optical emission lines which were observed to have a 
full width at half maximum of 200 to 2000 km s _1 . 

NGC 7213 hosts a blac k-hole with a mass of about 
10 s M (|Woo fc Urrvll2002ri. Its bolome tric luminosity of 
Lboi = 9xl0 42 erg s -1 ([starling et al.l2005t ) suggests a rather 
low accretion rate of 0.07 per cent of the Eddi ngton lumi- 
nosity (L E dd)- As noted bv lLobban et alj |2010r i. this value 
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is intermediate between those usually found in type I Seyfert 
galaxies and LINER's and it is significantly less than the pre- 
dicted 2 per cent LEdd 'critical' rate at which the 'soft state' 
transition appear s in black hole X-ray binaries (BHXRBs) 
IjMaccaronel I2003T ). Finally, based on its radio properties, 
NGC 7213 belongs to a rare class of extragalactic sources 
lying between the radi o loud and rad i o quie t AGN (it is one 
of the 20 sources of the lRov fc Norrisl (119971 ) sample). In this 
sense, it can be considered as an extragalactic analogue of 
the Galactic 'hard state' sources. 

The X-ray spectral behaviour NGC 7213 is also pecu- 
liar. IStarling et al] (|2005h using the reflection grating spec- 
trometer, on board XMM-Newton, found several emission 
features with no absorption lines; contrary to what is usu- 
ally observed in type I Seyfert galaxies. The absence of a 
Compton reflection component from either neutral or ion- 
ized material together with the lack of a relativistic Fe 
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Kq line suggests that the inner, optically thick acc retion 
disc in the source may be absent l|Lobban et al.l l2010), and 
replaced by an advec tion-dominated accretion flow (ADAF) 
|Starling et al.ll2005l ). This possibility further supports the 
'hard state' i nterpretation of t he source. 

Recently. iBell et al.l (|201 lh used radio and X-ray obser- 
vations of NGC 7213, from a long-term observing campaign 
consisting of several years, to search for correlated X-ray 
and radio variations that might originate in the so-called 
'fundamental plane of bla c k hole accretion' ((Merloni et al.l 
120031 : iKording et al.l I2OO6I ). IBell et all (|201lh showed that 
the average radio and X-ray luminosities fitted well with 
the global fundamental plane, consistent with a 'hard state' 
identification of this AGN. However, the X-ray versus ra- 
dio correlation within the monitoring period was weak, with 
significant intrinsic scatter away from the plane. 

Numerous studies in the past have shown that the X- 
ray photon index, V, correlates with the accretion rate, de- 
fined as £ = Lx/I/Edd (where Lx is usually the 2-10 keV 
X-r ay luminosity) in both AGN and BHXRBs. For exam- 
ple, IShemmer et all |2006l ) showed that F and £ are posi- 
tively correkited^jjsi^ a sam ple of thirty quasars, 
while [Sobolewska fc Papadakisl (|2009h found that the same 
positive correlation holds when one studies the spectral vari- 
ations of individual, luminous Seyfert galaxies. However, the 
positiv e F— £ corre l ation may not hold in less luminous AGN. 
In fact iGu fc Caq |2009l ) found an anti-correlation between 
r a nd £ in a sam p le of 55 low-luminosity AGN (LLAGN), 
and lYounes et all (|201ll ) reached a similar conclusion for a 
sam ple of 13 o p tically selected LINER's. 

IWu fc Gul (|2008i ) performed a detailed spectral study 
for six BHXRBs and found that F anti-correlates with £ be- 
low a 'critical' value of log(£ cr it) = —2.1 ± 0.2. At higher 
accretion rates, F and £ are positively c orrelated. Similar 
conclu sions were also reached recently bv lSob olcwsk a et al.l 
(2011). The r - £ relation in AGN may also change from 
negative to positive at a sim ilar ' critical' value, as it wa s 
suggested bv lWu fc Gul (|200Sl ). and lConstantin et al l J2009I ). 
This analogous behaviour between AGN and BHXRBs rein- 
forces the interpretation that LLAGN are analogues of the 
'hard state' BHXRBs while luminous Seyferts and quasars 
are the equivalent of the 'soft state' BHXRBs. On a first 
look, during the 'hard state', BHXRBs appear to have a 
constant hardness-ratio, in the commonly used hardness- 
intensity diagrams (e.g. iBelloni et all [2005). Nevertheless, 
detailed analysis of the lowest luminosity regime of this 
state, known also as 'low-hard' state, exhibits clearly an in- 
crease of the hardness-ratio for increased fl uxes, correspond - 
ing to an anti-correlation between F and £ (|Heil et al.ll2012h . 

The negative and the positive correlations between F 
and £, occurring below and above £ C rit, respectively, may 
be indicative for some 'switch' in the emission mechanism 
as the source's accretion rate increases above £ C rit- In fact, 
it has been suggested fWu fc Gul 120081 ; IConstantin et al.l 
l2009l ; lYounes et al.ll201ll ) that this transition could indicate 
the p assage from an ADAF (iNaravan fc Yi|[ l994; Esin et all 



Il997l ) to standard disc (but see Sobolewska et al 



(2011J) for 



alternative explanations as well). 

The positive correlation relation between F ver- 
sus £ in AGN has been firmly established statistically 
by either considerin g samples of numerous sources (e.g. 
Sh emmer et al.ll2006l ) or large data sets for a few individ- 



ual sources (e.g. llVLHardy et all 1 19991 ; lLamer et ail 120031 ; 
ISobolewska fc Papadakisl 120091 ) . However, up to now the F 
- £ anti-correlation has been established using on ly short X- 
ray observations of numerou s different sources jGu fc Cad 
120091 ; IConstantin et alj 120091 ; lYounes et ail l201ll ). while~it 
has never been observed in a single source. 

In this work, we use long-term, monitoring RXTE 
data, and present, for the first time, conclusive evidence for 
such an anti-correlation for an individual LLAGN, namely 
NGC 7213. We also use archival and proprietary data to con- 
struct its average spectral energy distribution (SED), and 
compare it with the average SED of luminous AGN and 
LINER's. In Section[2]we refer to the observations and data 
reduction procedures and in Section [3T] we present our re- 
sults consisting of the F versus flux relation and the spec- 
tral energy distribution of the source (SED) of the source. 
Finally, a discussion together with a summary can be found 
in Section 0] The cosmological parameters used through- 
out this paper are: Ho = 70 km s _1 Mpc -1 , f2A=0.73 and 
fi m =0.27, yielding a luminosity distance to NGC 7213 of 
22.12 Mpc (for a corrected redshift, z C orr.3K = 0.005145 to 
the reference frame defined by the 3 K cosmic microwave 
background radiation). This value of the luminosity distance 
appears to be fully co nsistent with the redshift-independent 
Tully-Fisher distance l|Tullvll 19881 ). 



2 OBSERVATIONS AND DATA-REDUCTION 
2.1 X-ray data 

NGC 7213 has been observed regularly by the Rossi X-ray 
Timing Explorer (RXTE) (proposal numbers: 92119, 93139, 
94139 and 94342) from 2006 March 3, 02:57:41 (UTC) to 
2009 December 30, 00:18:37 (UTC) (on-time: 800506 s). This 
is the complete set of X-ray observations of NG C 7213, ob- 
tained by RXTE, extending the X-ray data of IBell et ail 
|201ll ) by an additional of 100 days. In order to form the 
most homogeneous long-term data set, we use data obtained 
only by the proportional counter array 2 (PCU 2), one of the 
five Xenon pro portional counters forming the proportional 
counter array (|jahoda et al.l ll996). PCU0 and 1 have lost 
their propane layer, resulting an increased background rate 
as well as a different detector gain, and PCU 3 and 4 are 
usually switched off due to discharge problems. 

The Standard-2 (Std2: data with accumulation rate of 
16 sec) PCU 2 data (8 22 files in total) a re extracted and pro- 
cessed using f tools (|Blackburnlll995l ) included in HEAsof t 
(ver. 6.11.1), following the standard reduction procedures, 
provided by the RXTE Guest Observer Facility (GOF). For 
each observation we create a filter file, assembling the scien- 
tifically important parameters, using the script xtefilt and 
then by using the f tool fmerge we merge them all in a sin- 
gle 'master' filter file. Then, based on this file, we extract 
the useful observing time periods, known as good time inter- 
vals (GTIs), using the f tool maketime having the following 
observational constraints: elevation angle greater than 10°, 
a pointing offset of less than 0.02°, electron contamination 
less than 0.1, time since south Atlantic anomaly between 
zero and thirty minutes and time since PCU breakdowns less 
than -150 s or greater than 600 s. Finally, for each Std2 ob- 
servation a synthetic background model file is created with 
the ftool PCABACKEST using the 'faint background model'. 
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For the production of the background-subtracted light 
curves, in the 2-4, 5-10 and 2-10 keV energy bands, 
we use the ftool SAEXTRCT to produce the 'source-plus- 
background' and 'background' light curves for each energy 
band, respectively. We select events obtained during the 
GTIs registered only from the top xenon layer (X1L, X1R) 
of PCU 2, optimizing in this way the signal to noise ratio. 
Then the corrected background-subtracted light curves are 
produced by using the ftool lcmath in units of count-rate. 
In order to produce a total long-term X-ray light curve in 
the 2-10 keV energy range, in units of flux i.e. erg s _1 
cm -2 , we extract for each Std2 observation, and the cor- 
responding synthetic background file, a spectrum together 
with the PCU 2's response using again the ftool SAEXTRCT 
and the script PCARSP, respectively. Again, we select events 
obtained during the GTIs which are registered only from 
the top xenon layer of PCU 2. The resulting spectra are then 
grouped so that each spectral bin contains at least 30 counts, 
using the ftool task GRPPHA. 

Usin g the X-ray fitting package XSPEC (ver. 12.7.0) 
(|Arnaudl 11996) we fit to each X-ray spectrum a power 
law model assuming photoelectric absorption (XSPEC model 
wabs), of a fixed interstellar column-density of JVh = 
1.1 x 10 20 cm" 2 (e stimated using the FTOOL nH, after 
iKalberla et al.l 120051 ) . Finally, after fixing the fitting pa- 
rameters (i.e. photon index and normalisation) to their the 
best-fitting values, we obtain the flux values, in the 2-10 
keV energy range using the convolution model cf lux. The 
errors in the spectral indices and fluxes indicate their 68.3 
per cent confidence range, corresponding to a A\ 2 of 1, un- 
less otherwise stated. 



2.2 Spectral energy distribution data 

In order to construct the SED of NGC 7213, we anal- 
yse long term monitoring radio and optical data. Below 
we describe the data analysis procedures for these obser- 
vations. We also use the p u blished near-infrared , nuclear 
fluxes of iHonig et~ai1 ( |201Ch ; lAsmus et al.l (|201lh at 12.27 
and 11.25 and 10.49 /im (these fluxes are listed in Table 
[T]with an asterisk). Finally, we also considered the results 
from a 130 ksec XMM-Newton observation obtained during 
55146-55148 MJD (obs: ID-0605800301, Emmanoulopoulos 
et al. in prep.), the 4.5 years average Swi ft-BAT spectrum 
(as r eported in entry 1185 of table 2 in ICusumano et all 
|2010T), and the recen t 0.1-100 GeV Fermi-L AT upper limit 
l Lobban et"ai]|20ld ). 

2.2.1 Radio data 

We searched the Australian Telescope Online Archive 
(ATOA) for Australian Telescope Compact Array (ATCA) 
observations of NGC7213's location, selecting a range that 
covered as many frequencies as possible, opting for the 
longest duration observation (s) in each case (sometimes 
yielding multiple useful files, as in the case of 5-10 GHz). 
While attempts were made to process all designated files, 
the highest available frequencies, greater than 20 GHz, did 
not yield useful images. The final data set used in this work 
was comprised of 21 ATCA observations, from four projects; 
C782, C1803, C1532 and C1392. 



Table 1. The radio, near-infrared and optical mean flux values 
of NGC 7213. 



Energy band Flux 
mjy 



1.344 GHz 


121.3 ± 2.2 


1.384 GHz 


112.4 ±9.3 


2.386 GHz 


114.9 ±3.6 


2.496 GHz 


98.1 ±4.1 


4.8 GHz 


135.8 ±8.3 


8.64 GHz 


150.6 ± 24.3 


17 GHz 


140.4 ±3.9 


19 GHz 


128.5 ±4.1 


12.27 Atm* 


235.8 ± 128.4 


11.25 Atm* 


232.9 ± 144.5 


10.49 ^m* 


239.1 ± 22 


5500 A 


0.49 ±0.18 


4400 A 


0.73 ±0.15 



* the n ear-IR measu r ement s taken from [Honig et af, 
lAsmus et al.1 j201 lh as explained in the text. 

The primary calibrator for all these data sets was 
PKS 1939-6342 (PKS B 1934-638). C782 includes observa- 
tions from 1999 March 03 and frequencies of 1344 MHz 
(0.3 h on source, average beam size (ABS): 64.7" x 7.0"), 
1384 MHz (8.5 h on source, ABS: 21.6" x 7.2") and 2496 
MHz (8.8 h on source, ABS: 7.6" x 4.9" ), with the anten- 
nae in configuration 6C (minimum baseline length of 153 
m, maximum of 6 km). The secondary calibrator used was 
PJS J2214-3835 (PKS B2211-388A). C1803 included obser- 
vations from 2008 April 21 and 23 at 1384 MHz and 2386 
MHz, with 0.15 and 0.1 h on source respectively. The an- 
tenna configuration was 6A (minimum baseline length of 337 
m, maximum of 5939 m) and the secondary calibrator used 
was PKS J2235-4835 (PKS B2232-488). C1532 contributed 
all our 4800 MHz (ABS: 17.3" x 1.8") and 8640 MHz (ABS: 
9.64" x 1") data, with observations covering the dates of 2008 
January 05, 24, April 12, November 04 and 20 with on source 
durations of 1.48, 1.7, 1.5, 2.29 and 1.49 h, respectively. 
The secondary calibrator used in all cases was PKS J2218- 
5038 (PKS B2215-508) and antennae were always in config- 
uration 6A. Finally, C1392 provided our highest frequency 
ATCA data, at 17 GHz (ABS: 35.5" x 27.55") and 19 GHz 
(ABS: 29.15" x 22.8"). Two observations are used, from 2009 
September 29 and 30 with 0.02 and 0.03 h on source re- 
spectively. The secondary calibrator was PKS J2248-3235 
(PKS B2245-328) and the antenna configuration was H75 
(minimum baseline of 31m, maximum of 4408m). 

Once imaging was complete, point-source fitting is used 
to measure flux density at NGC 7213's position RA(1950) = 
22 h 09 m 16*26, Dec. (1950) = -47° 09' 59"95. All data and 
image processing is carried out in the radio interfe rometry 
data reduction package MIRIAD (|Sault et alj|l995h . 



2.2.2 Optical data 

We used the ANDICAM instrument mounted on the 1.3 m 
telescope SMARTS in Chile to observe NGC 7213 in opti- 
cal B (4400 A) and V (5500 A) bands. Observations were 
performed between August 2006 and August 2011, every 4 
days and during each observations two 60 s exposure frames 
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in the B band and two 30 s exposures in the V band were 
obtained when the target was visible. After discarding the 
'bad seeing' nights, a total of 266 epochs are accumulated in 
each band. The field of view of ANDICAM is 6' x 6' covering 
both the entire host galaxy and a few bright stars, and the 
pixel size is 0.37" per pixel. The typical seeing during the 
campaign was about 1.5". The data reduction for both the 
B and V band images is the same. 

Since NGC7213 host galaxy is a face-on spiral with a 
bright nucleus compared to the AGN, galaxy subtraction has 
to be performed to determine the pure AGN flux. To model 
correctly the large scale struct ure we use the d ata analy- 
sis algorithm GALFIT (ver. 3.0) (jPeng et al]|2010h . fitting a 
Nuker function. After the fitting, a bright ring 10" from the 
nucleus is left together with some spiral structure inward of 
the ring and a resolved core. The core is subtracted using 
a Moffat function leaving a residual nuclear peak consis- 
tent with the stellar point-spread function (PSF), which we 
identified with the AGN. Aperture photometry is performed 
with the task da ophot, of the IRAF (ver. 2.15) software sys- 
tem l|Todvlll993h . on the stars of the original image and on 
the modelled nuclear PSF to determine its relative bright- 
ness. In order to estimate the flux, the reference stars in turn 
are calibrated taking their average aperture magnitudes over 
all reported photometric nights when observations were car- 
ried out, corrected for airmass extinction. The average nu- 
clear fluxes that are obtained in this way are 5 x 10~ 12 erg 
s" 1 cm" 2 at 4400 A and 3 x 10" 12 erg s" 1 cm" 2 at 5500 
A, while the fractional root mean square variability ampli- 
tude (including observational noise and intrinsic variability) 
is around 1 erg s" 1 cm" 2 for both bands, respectively. 



3 RESULTS FOR NGC 7213 

3.1 The X-ray flux — photon index relation in 
NGC 7213 

The upper panel of Fig. [T] shows the 2-10 keV, X-ray light 
curve of NGC 7213. One can observe variations on time 
scales of days and weeks, superimposed on a flux decreasing 
trend, from the start to the end of the monitoring campaign. 
Each measurement on this plot corresponds to one RXTE 
pointing, having a mean duration of 1.06 ±0.25 ks (the error 
estimate corresponds to the standard deviation of the ex- 
posure times of individual observations). On average RXTE 
observed the source every ~ 2.3 days, except from the period 
between 54982-55077 MJD (indicated by the double-headed 
arrow in the upper panel of Fig.[U), where observations were 
performed around every half a day. Observations with large 
uncertainties e.g. around (53797+450) MJD, correspond to 
exposure times less than 0.9 ks. 

The lower panel of Fig. [T] shows the evolution of the 
F as a function of time. In order to investigate any possi- 
ble relationship between flux and F, we group both data 
sets in bins of 45 consecutive observations (each bin is then 
~ 100 days long). For the intense monitoring period (1185- 
1280 MJD-53797) we apply the same sampling scheme by 
selecting only the observations that are separated around 
2.3 days and then bin them in bins of 45 consecutive ob- 
servations. In this way we ensure that the resulting X- 
ray flux and photon index range are represented by equal 




200 400 600 800 1000 1200 1400 
Time (MJD-53797) 

Figure 1. The RXTE results of NGC 7213. (Top panel) The long 
term light curve in the 2—10 keV energy range. The double-headed 
arrow indicates a period intense monitoring between 54982-55077 
MJD. (Lower panel) The photon index of each pointing observa- 
tion as a function of time. 

number of observations, avoiding the possibility that the 
larger number of observations in the intense monitoring pe- 
riod could drive the resulting flux/photon index relation. 
We estimate a weighted me an and a weighted error (e.g. 
iBevington fc Robinson! Fl992) for both flux and F in each 
bin. Fig. [2] shows the resulting average X-ray flux and V 
values, plotted against each other, unveiling a clear 'harder 
when brighter' behaviour. 

In order to quantify this anti-correlation w e first com- 
pute the Kendall's r rank correlation coefficient (|Press et al.l 
1992), using the data plotted in Fig. yielding a r = —0.81 
with a null hypothesis probability of 2.6 x 10" 5 . This re- 
sults shows that the anti-correlation between spectral slope 
and flux is highly significant. Then, we fit to the data a lin- 
ear model, (y = ax + 0), consider ing the uncertainties in 
both coordinates (|Press et al.l 1992). The best-fitting model 
yields a slope of a — —0.063 ± 0.018 and an intercept of 
P = 1.99 ± 0.04 with a X 2 merit function of 2.90 for 13 de- 
grees of freedom (d.o.f.) having a null hypothesis probabil- 
ity of 1.7 x 10~ 3 . Fig. |3]shows the confidence contour (solid 
black line) used to estimate the one-standard deviation un- 
certainties on the best-fitting parameters corresponding to 
an ellipsoidal region with a \ 2 of 3.90 (A^ 2 = 1 from the 
minimum for 1 d.o.f.). In the same plot, we show also the 
contours for confidences of 95 and 99 per cent as well as the 
68.3 per cent joint confidence contour for both the slope and 
the intercept with the latter corresponding to a A\ 2 = 2.30 
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Figure 2. The average photon index versus the average X-ray 
flux of NGC 7213 in bins of 45 consecutive observations. The data 
suggest a 'harder when brighter' relation. The solid line indicates 
the best-fit linear model to the data taking into account the un- 
certainties in both coordinates. 



from the minimum for 2 d.o.f. A simple linear regression 
model, taking into account only the photon index uncer- 
tainties, yields equivalent results: a = —0.063 ± 0.008 and 
P = 1.99 ± 0.02 (x 2 =2.90 for 13 d.o.f). The resulting X 2 , 
from both methods, are relatively small, indicating that the 
estimated errors on the average T within each bin maybe 
slightly overestimated. We therefore repeat the fit and this 
time we fit to the data a linear model following th e 'ordinary 
least-s quares regression of Y on X' routine of llsobe et al.l 
(1990) which does not take into account the error on the 
data. The best-fit results from these routines are consistent, 
within the quoted errors, with the previous results yielding: 
a = -0.072 ± 0.008 and j3 = 2.00 ± 0.02. 

We have performed several sanity checks in order to 
test the sensitivity of our results to the binning scheme. 
We have considered alternative binning of 20,60 and 100 
consecutive observations as well as considering all the data 
from the intense monitoring period. Also we have per formed 
data-resampling using jackknifing (Sh ao fe Tu||l995T ) by se- 
lecting randomly sub-samples from the data set and binning 
them in the flux/photon index plane. All the results are ex- 
tremely consistent with each other. Finally, in order to en- 
sure that the overall anti-correlation trend is not induced 
by a ny sort of statistical depen dence between the flux and F 
(e.g. IVaughan fc Edelsonll200ll ), we measured the eccentric- 
ity and the orientation of the \ 2 contour plots of the fits, 
used to derived the uncertainties of the fluxes and slopes. 
The 68 per cent ellipsoids have an average eccentricity of 

0. 91 ± 0.03 and they are tilted on average by an angle of 
(79 ±2)° counter-clockwise from the horizontal (flux) axis 
favouring, if at all, trends moving the opposite direction 

1. e. a positive correlation between flux and V. We therefore 
conclude that for NGC 7213, contrary to what is observed 
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Figure 3. The confidence contours for the best fitting model of 
Fig. [2] taking into account the uncertainties in both coordinates. 
The black thin lines correspond to the 99, 95 and 68.3 confidence 
contours and represent the ellipsoids having from the minimum 
a Ax 2 of 6.63, 3.84 and 1, respectively, for 1 d.o.f. The grey 
thick line corresponds to the the 68.3 per cent joint confidence 
contour for the slope and the intercept having from the minimum 
a Ax 2 of 2.30, for 2 d.o.f. The horizontal and vertical dashed grey 
lines correspond to the tangents of the 68 per cent confidence 
region yielding the 68 per cent uncertainty ranges of the slope 
and intercept respectively. 



in other luminous Seyfert galaxies, the X-ray photon index 
anti-correlates with the X-ray flux. 



3.2 Hardness Ratio Analysis 

A completely model independent way to check for the valid- 
ity of the above-mentioned anti-correlation behaviour, is by 
estimating the hardness ratio from the X-ray light curves. 
After binning the light curves in bins of 50 consecutive ob- 
servations, we estimate the hardness ratio (5-10 keV)/(2-4 
keV) versus the overall count-rate in 2-10 keV. In Fig. [4]wc 
plot the corresponding estimates, showing a clear increasing 
trend which implies that the X-ray behaviour of the source 
becomes 'harder' when the source gets brighter. 

We can now fit to the hardness ratio data a linear model, 
(y — ax + f3), taking into account the errors in both coordi- 
nates (as in Section [3T}. The best fit model has a slope of 
a = 0.11 ± 0.01 and an intercept /3 = 1.81 ± 0.02 yielding a 
X 2 of 4.31 for 15 d.o.f. having a null hypothesis probability 
of 3.4 x 10 -3 . Since the 99 per cent confidence interval for 
the slope is (0.07,0.14) the null hypothesis can be rejected at 
1 per cent significance level. Finally, the value of t-statistic 
that we get from the data is 10. iQ, corresponding to a prob- 
ability of getting such a value from chance alone equal to 



1 The value of t-statistic is ti5,o.o05 = 2.95 
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Figure 4. Hardness ratio plot of the count-rate in (5-10 keV) / (2- 
4 keV) versus 2—10 keV. The bin size corresponds to 50 consecu- 
tive observations. The solid line indicates the best-fit linear model 
to the data taking into account the uncertainties in both coordi- 
nates. The positive slope implies a 'harder when brighter' X-ray 
spectral behaviour. 



4.6 x 10 . Therefore, we can robustly conclude that the 
best-fit slope is significantly different from zero. 



3.3 The nuclear broad-band Spectral Energy 
Distribution 

The average spectral energy distribution of NGC 7213 (in 
vL v representation) is shown in Fig. [S] (both panels) with 
black symbols corresponding to flux estimates (listed in Ta- 
ble [1} derived in this work and and grey symbols to archival 
data. The error bars correspond to the standard deviation 
around the mean flux, at a given frequency, whenever multi- 
ple observations from different epochs are available. In this 
way, the error bars are indicative of the actual variability of 
NGC 7213 (the actual uncertainty of single pointing obser- 
vations is much smaller than the symbol size). The average 
flux from the ensemble of RXTE observations in the 2-10 
keV energy range is shown by symbol 1, and corresponds to 
a power law having normalization of (6.51 ± 3.73) x 10~ 3 
photons s _1 cm -2 keV -1 (at 1 keV) and a photon index of 
1.87T0.23. The bowtie 2 indicates the 0.1-10 keV spectrum 
from the XMM-Newton observation (Emmanoulopoulos et 
al. in prep.) which is broadly consistent with a power law 
with a normalization 2.99^qq3 x 10 -3 photons 
keV -1 (at 1 keV) and a photon index 1.6 
3 indicates the average flux estimate registered by the Swift- 
BAT instrument, between 15-150 keV, corresponding to a 



cm 

1°;°*. The bowtie 



power law with normalization^ 5.4l 3 'J x 10 3 photons s 1 
cm -2 keV -1 (at 1 keV) and a photon index of 1.82lto'ia- 
Finally, the 0.1-100 GeV the Fermi-LAT upper limit, as- 
suming a photon index of 1.75, is 3 x 10~ 9 photons s _1 
cm" 2 , is indicated by the arrow at 0.25 GeV (mean energy 
of the Fermi-LAT energy band). 

The mean SED of NGC 7213 is shown in Fig. [S] This 
is true even for the RXTE and Swift-BAT data, since at 
these flux levels the contribution of any host galaxy emis- 
sion should be minimal, especially in the Swtft-HAT band, 
as well as for the optical data, where we have carefully dis- 
card the host galaxy contribution from our measurements. 
Also note that, the optical, RXTE and Suii/i-Bat data plot- 
ted on this figure are indicative of the source flux averaged 
over many years, which are largely overlapping. On the other 
hand, the beam size of the radio observations is rather large, 
due the short exposure times (Section 12. 2) . As a result, the 
observed radio fluxes we have used may be contaminated by 
radio emission from a substantial part of the host galaxy 
itself. However, our radio flux estimates, at least for the 
4.8 and 8.6 GH z bands, ar e iden tical to the mean flux val- 
ues reported bv lBell et al.l l|201ll ). who used a beam size of 
0.5" for a 12 h integration time. In addition, the observed 
variability at all radio bands indicate that most of the emis- 
sion should originate from the nucleus itself. Detailed X-ray 
studies (e.g. Emmanoulopoulos et al. in prep.. [Cobban et al.l 
2010) indicate that the intrinsic absorption towards the nu- 
clear source in NGC 7213 is minimal in addition to the fact 
that the galaxy of NGC 7213 is face-on. 

Consequently, the SED in Fig. [5] should be representa- 
tive of the mean-intrinsic nuclear SED of NGC 7213. The 
blue solid and dashed-line in the top panel of Fig. [S] indi- 
cates the average SED of radio-quiet and radio-loud quasars 
l|Elvis et al.lll994T ). respectively, and the blue open-circles in 
the bottom pane l of the same figure i ndicate the average 
SED of LINERs l|Eracleous et alj|201(f ). Their relative po- 
sition is set by minimizing the distance of their logarithmic 
ordinates from those of NGC7213's, at the same frequen- 
cies, weighted by their squared errors. Since the distance 
minimization procedure is done in the logarithmic plane, in 
this way we estimate effectively an optimum normalization 
for each average-SED. Our results suggest that it is the LIN- 
ERs mean SED which fits best, within the NGC 7213's SED. 
We therefore conclude that the nuclear continuum emission 
of the source is indeed representative of the mean SED of 
nearby LINER's. 

Despite the fact that the (archival) near-IR flux mea- 
surements of NGC 7213 (as shown in the bottom panel of 
Fig. [5} are representative of emission from a central region 
which is less than 0.35" in size, it is not clear whether they 
correspond to emission form the nuclear source itself. For 
example, if we indeed observe the active nucleus directly, 
without any significant absorption, the detected near-IR 
emission cannot be due to optical/UV nuclear emission be- 
ing absorbed by the putative dusty, obscuring torus in this 
galaxy and re-emitted in the near-IR, because the average 
optical luminosity is significantly lower than the observed 



The flux value is consistent with the one derived from the 



maps 



of th e 4 th IBIS/ISGRI soft 7-ray survey catalogue llBird et al.l 
l201Ct> . 
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Figure 5. The SED of NGC 7213. For both panels, black sym- 
bols correspond to flux density estimates derived in this work 
and grey symbols correspond to archival data. The enumerated 
symbols, 1,2 and 3, correspond to the average RXTE spectrum 
(depicting the long-term variations), XMM-Newton and Swift- 
BAT spectrum, respectively. The grey arrow, corresponds to the 
Fermi-LAT upper-limit, in the 0.1-100 GeV energy range. (Top 
panel) The SED of NGC 7213 together with the average SEDs of 
radio-quiet and radio-loud quasars shown w ith the blue solid - 
and dashed-lines, respectively (taken from lElvis et al.l Il994l) . 
(Bottom-panel) The SED of NGC 7213 together together with 
the average SEDs LINERs shown with blue open-circles (taken 
from lEracleous et al.ll20ld ) . 



near-IR luminosity. On the other ha nd, NGC 7213 shows 
evidence of a starburst driven wind ( Bian chi et al.l 120081 ) 
and hosts a star-forming ring few kpc awa y from the nu- 
cleus . A clumpy torus with toroidal shape (|Nenkova et al.l 
2008), having a small angular width parameter can then 
explain the observed infrared luminosity as being due to 
reprocessed radiation from star-forming regions, while the 
optical nuclear emission could escape unabsorbed through a 
torus-hole. Note that LINERs are thoug ht to have rela tively 
hot yet normal main-sequence O stars (|Shields!ll992l ) able 
to heat the dust around them, something that strengthens 
even more the LINER interpretation of NGC 7213 SED. 

The SED plotted in Fig. [S] is based on more observa- 
tions of the nuclear source flux, and is sprea d over a larger 
freque ncy range, than the SED presented by Istarling et al.l 
(2005). We can therefore use it to derive a more accurate 
estimate of the nuclear bolometric luminosity. After inter- 
polating the SED linearly in logarithmic space, we integrate 
it between 1.344 GHz and 3.63 x 10 19 Hz (i.e. 150 keV) and 
we derive Lboi = 1.7 x 10 43 ergs s _1 , yielding an accretion 
rate of 0.14 per cent of Eddington limit. 



4 DISCUSSION 

We have analysed the long-term RXTE observations of 
NGC 7213 and we found that this low luminosity source 
exhibits a clear 'harder when brighter' X-ray behaviour 
(Fig. [2]). This is the first time that such a spectral variability 
behaviour is reported for either a 'low' or 'high' luminos- 
ity AGN, and is contrary to what is observed in luminous 
Seyferts and quasars. We also constructed the average, nu- 
clear SED of this source, using archival and proprietary data 
that we reduced, together with other measurements from 
the literature and we found that its shape closely resembles 
the mean SED of nearby LINERs. Finally, we provided a 
new measurement for the nuclear bolometric luminosity of 
NGC 7213: Lboi = 1.7 x 10 43 ergs s _1 , yielding a rather low 
accretion rate of 0.0014 (i.e. 0.14 per cent) of the Eddington 
limit. 

The average, broad band nuclear SED of NGC 7213 re- 
sembles that of LINERs (Fig.[5J, bottom panel). In the aver- 
age SEDs of radio- loud and the radio-quiet AGN, the optical 
flux is quite higher than the X-ray flux. This is definitely not 
the case for this source, since the optical flux lies well be- 
low that of the X-rays. This is a very robust result, as both 
the optical and the X-ray data are based on long, monitor- 
ing observations, i.e. they show the ave rage behaviour of the 
source over many years. Furthermore, IWu et all j 19831 ) no- 
ticed that the UV excess for NGC 7213 should be extremely 
week or absent, while both the recent XMM-Newton obser- 
vations (Emmanoulo poulos et al, in prep .), as well as the 
Suzaku observations (jLobban et all |2010T ) show no indica- 
tion of soft X-ray excess, which is typical in most luminous 
AGN. All these results strongly suggest that the 'big-blue 
bump' is missing in this source. This can be expected in a 
scenario in which the inner part of the geometrically thin 
and optically thick disc is missing in NGC 7213. 

Our results suggest that the accretion rate of the source 
is significantly smaller (by a factor of 10) than the 'criti- 
cal' rate at which accreting BHXRBs move from the 'hard' 
to the 'soft state'. Therefore, NGC 7213 could be the 'hard 
state' analogue of BHXRBs, and in fact, its 'harder when 
brighter' behaviour strongly supports this hypothesis. Al- 
though the global relationship between F and £ in BHXRBs 
is well-established by c omparing measurements from single - 
epoch observations (e.g. lWu fc Gull2008l ; lYo"unes et al.ll201ll ) 
the correlation of V with £ on short time scales (within an 
observation, i.e. down to minutes or even seco nds) is less 
well-determined. However. I Axelsson et all ((2008;) showed for 
the BHXRB CygX-1 that the hardness-flux anti-correlation, 
seen in variations within an observation in brighter 'hard' 
states, becomes a positive correlation in the faintest 'hard' 
states. This change in behaviour is consistent with a positive 
r-£ correlation at higher £ changing to an anti-correlation 
at lower £, i.e. the variations within an observation follow 
the same trends as the global F-£ relationship. 

Theoretically a considerable progress has been accom- 
plished in this field with a better understanding of the com- 
plexity of the BHXRBs' 'har d state' as well as LLAGN (see 
for a review, iNaravaril 120051) . A currently proposed model 
involves an accretion disc plus a hot accretion flow model, 
ADAF model, to explain the spectral dependence on accre- 
tion rate for BHXRBs, from quiescence up to the 'soft state'. 
For intermediate accretion rates, as the accretion rate in- 
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creases the Compton parameter in the hot accretion flow 
increases as well, producing a Comptonisation X-ray spec- 
trum with a 'harder' power law slope. This regime is iden- 
tified with the 'hard state', at mass accretion rates of 1-8 
per cent of the Eddington rate, since the accretion process 
there is inefficient corresponding to lower fractions of LEdd, 
consistent with our observations of NGC 7213. As the ac- 
cretion rate increases further, the model predicts a differ- 
ent behaviour in which increasing luminosity corresponds to 
'softer' X-ray power law slopes. This regime is identified with 
the 'soft state' and this is the what is normally observed in 
higher luminosity AGN. 

Finally, another interesting feature of the source's SED 
is the fact that its X-ray emission above 20 keV is quite 
high, and implies a high-energy cut-off larger than 350 keV 
i|Lobban et al.ll2010h . NGC 7213 shows some week evidence 
of a radio jet structure at 8.4 GHz dBlank et all |2005| V If 
such a structure is indeed there and it is aligned towards 
the observers direction it could produce relativistically am- 
plified radio emission through synchrotron radiation, and en- 
hanced X-ray emission through inverse Compton radiation 
of either the synchrotron photons and/or optical photons 
from the star-burst environment of the host galaxy. This 
possibility could explain the fact that the radio emission of 
the source lies between that of radio-loud and radio-quiet 
AGN, similar to blazars. At the same time, the existence 
of a jet could also account for the 'harder when brighter' 
behaviour of NGC 72 13, something which is commonly ob- 
served in blazars (e.g. | Krawczvnski et al 1 120041 :151 iozzi et al.l 
2006; Zhang et al.ll2006h and can be explained in the frame- 
work of the synchrotron self Compton models. A jet model 
has also been proposed to explain the spectral evolution of 
the BHXRB XTEJ1550-564 as it moves from the fad ing 
state towards the X-ray 'hard state' (|Russell et alj|201(ih . 

Therefore, a jet pointing towards us, contributing sig- 
nificantly to the radio and X-ray emission of the source can 
not be ruled out, and in fact, for the case of ADAF mod- 
els, it can be created through the B landford-Znajek mech- 
anism (|Armitage fe Nataraianlll999l ). In the 'hard state' of 
BHXRBs the kinetic lumi nosity in the jet is be lieved to equal 
the radiation luminosity l|Fender et alj[201ol 'l. For the case 
of NGC 7213, assuming the existence of a jet, the previously 
derived accretion rate may therefore be lower by a factor 
of two, still remaining well below the typical transition be- 
tween the 'hard' and the 'soft' state in BHXRBs. In the 
future, models combining a j et component wit h a geomet- 
rically thin disk and ADAF (|Yuan et all 120051 ') need to be 
tested for the case of NGC 7213. 
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